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Abstract 


We define the generalized hypergeometric polynomial of degree N as follows: 


N 


P N (a 1> ...,a p ;P 1 ,...,P q ;z) = ^ 


m —0 


KLhiL-'KL z N ~ m 


m[ (Pl)m • ' ' Pd 


9 / m 


= z N p+1 F q (-N,a p ; (3 1 ,..., p q , 1 /z) 


Here N is an arbitrary positive integer, p and q are arbitrary nonnegative integers, the p + q parame¬ 
ters a.j and j3 k are arbitrary (’’generic”, possibly complex) numbers, (a) m is the Pochhammer symbol and 
p+iF q (ao, cci, ...a p ; P 1 , ..., fi q ; z) is the generalized hypergeometric function. In this paper we obtain a set of 
N nonlinear algebraic equations satisfied by the N zeros of this polynomial. We moreover manufacture an 
N x N matrix L in terms of the 1 + p + q parameters N, ay, (3 k characterizing this polynomial, and of its 
N zeros ( n , and we show that it features the N eigenvalues A m = m nLiHh + i-m) , to = 1,..., IV . 
These N eigenvalues depend only on the q parameters (d k , implying that the N x N matrix L is isospectral for 
variations of the p parameters ay; and they clearly are integer (or rational) numbers if the q parameters p k 
are themselves integer (or rational) numbers: a nontrivial Diophantine property. 

Keywords: hypergeometric polynomials, diophantine properties, Jacobi polynomials, isospectral matrices, 
special functions. 

MSC: 33C20; 11C08; 11C20. 


1 Introduction 

The investigation of the properties of the zeros of polynomials has a long history going back for several 
centuries, yet new approaches and findings have also emerged in relatively recent times. A class of such findings 
extended the results pioneered by G. Szego (see in particular Section 6.7 of CD, by identifying additional sets 
of nonlinear algebraic relations satisfied by the zeros of the classical polynomials and, more generally, of 
polynomials belonging to the Askey scheme, as well as N x N matrices, constructed with the N zeros of 
these polynomials (of degree N), whose eigenvalues could be explicitly identified and in many cases feature 
Diophantine properties: see for instance m 0 a 0 m 0 m- The present paper reports analogous results 
for generalized hypergeometric polynomials. These findings are displayed in the following Section 2, and proven 
in the subsequent Section 3; certain polynomial identities which are essential for obtaining and reporting 
these results—and are themselves remarkable—are proven and displayed in Appendix A. A terse Section 4 
(’’Outlook”) outlines possible future developments. 
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2 Results 


The generalized hypergeometric function p +iF q (ao, <+i, • ••, a p ; /3 lt ..., (3 q \ z ) is defined as follows (see for instance 

my- 


p+1 Fq (ao, «1, •**) /f^, ..., /5^, z) ^ ' 

J=0 


(«o)y («i)y • • • (a P )y ^ 

j! (Pi)j ■ ■ ■ (Pq), 


(la) 


Above and throughout, the Pochhammer symbol (a) ■ is defined as follows: 


(a) 0 = 1 ; (a) - = a (a + 1) • • • (a + j - 1) = 


r (u + j) 
P(a) 


for j = 1,2,3,.... 


(lb) 


Clearly if one of the p +1 parameters ay is a negative integer, say ao = —N, and all the other p + q parameters 
ay and /3 k have generic (possibly complex) values, the series in the right-hand side of the definition (fTal) of 
the generalized hypergeometric function terminates at j = N (since (—IV) . = 0 for j = N + 1, N + 2,...) 

. Hereafter we call generalized hypergeometric polynomial the resulting polynomial (of degree N in z, and 
conveniently defined as follows, so that it is monic): 


N 

P N (ai,...,ot p ;f3 1 ,...,/3 q ;z) = ^ 

m —0 


HLK-KL z N ~ m 


i - (Pi) 




QJ m 


(2a) 


and we denote its N zeros as ( n , 

Pn (a 1 ,-,a p ;fi 1 ,...,/3 g ;( n ) = 0 , n=l,...,N. (2b) 

Hence the values of the N numbers ( n depend on the 1 +p + q parameters N , ay, f3 k . 

Notation 2.1. Above and hereafter N is an (arbitrarily assigned) positive integer, p and q are two (arbitrarily 
assigned) nonnegative integers, and indices such as n, to, £ (but not necessarily j, k ) run over the N integers 
from 1 to TV (unless otherwise indicated). The N zeros £ n are of course defined up to permutations. In the 
following we always assume the same assignment to be made for the correlation of the values of the N zeros 
of a polynomial with the values of the index n labeling them. And below underlined lower-case letters denote 
A^-vectors (hence, for instance, ( = (Cd Cjv))i an d underlined upper-case letters denote N x N matrices 
(hence for instance the matrix L has the N 2 elements L nm ). Finally: we always adopt the standard convention 
according to which a sum containing no terms vanishes, and a product containing no terms equals unity: for 
instance, £j,fc=i u yfc = 0, H; A - I:y/A - = L D 

The first result of this paper consists of the following 

Proposition 2.1. The (set of) N zeros Q n of the generalized hypergeometric polynomial Pn (a i, ..., a p ; /3 1; ..., /3 g ; z), 
see (P2al) and m, satisfy the following system of N algebraic equations: 


9+1 r 




= 0 


3=0 


n = 1,..., N . 


( 3 ) 


Here the q + 1 coefficients b k , respectively the p+1 coefficients ay, are defined in terms of the q parameters /3 fc 
respectively the p parameters ay so that 


9+1 


fc=l 


hence 


fc=l 


6 2 - - E 

3 =1 


E ( bk xk ) > 

(4a) 

fc=i 


-i) , 

(4b) 

(.Pk ~ 1) ) 

(4c) 


k= 1, kzfij 
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and so on up to 
respectively 

hence 


and so on up to 


<*=5 i 


ej=i-e^j L fc=1 ’ 

b q+ 1 = (-1) 9 ; 


n (/?* -1) 


V x , 

n k- - z) = e a o xJ 


7 =1 


3=0 


ao= n ( a j) . 

7=1 


v 

ai = ~j2 

k =1 

1 V 

^ = \ E 

l,k=l;l^k 


I! («j) 

j'=i; 


n («i) 

7=1; jW.fe 


a p = (-1 ) P ■ 


As for the functions fn^ (C) of the A” zeros £ m , 


they are defined recursively as follows: 


E' +1) (c) 


N 

-fS ) «J+ E 

£=1; 


c„/j i) (a+c//n o) (a 

Cn 0: 


J — 1,2,3,... , 


(4d) 

(4e) 

(5a) 

(5b) 

(5c) 

(5d) 

(5e) 

(6a) 


with 

fi 1] (C) = Cn - (6b) 

implying the expressions of fn ^ (C) with j = 1, 2,3,... reported in the Appendix, see (l59l) . 

And the functions g^} (£) of the N zeros f n are defined as follows: 

9 { n ] (C) = 1 , (7a) 

N 

9n ] (C) = E 

£=1; l^n 

implying 



and the expressions of g£* (£) with j = 1, 2, 3,... reported in the Appendix, see m - D 

Remark 2.1. The functions fn ^ (C) and gffl (C) are universal: they do not depend on the generalized 
hypergeometric polynomial under consideration. But of course their arguments do, being the N zeros f n of 
the polynomial Pn («i, ..., a p \ /3 1 ,..., /3 ? ; C) , see ©. □ 

The second (and main) result of this paper consist of the following 

Proposition 2.2. Let the (unordered) set of N numbers £„ denote the N zeros of the generalized hyperge¬ 
ometric polynomial Pn («i,..., a p \ (3 1 , ..., (3 q \ z ), see (j^al) and (I2bl) : and let the N x N matrix L (£) be defined 
componentwise as follows, in terms of these N zeros and the 1 +p + q parameters N, ctj, j3 k characterizing the 
generalized hypergeometric polynomial Pn (ai,..., a p ; ^ l5 ..., /3 ? ; z): 

9+1 P 

Lnm (C) = E [ bk fmm (C) 

fc=l 7=1 


E 


«7 9nL 


(C) 


( 8 ) 


W (C) + fP (C) 

Cn Cl 


3 = 1,2,... 


(7b) 
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where of course the coefficients b k and aj are defined as above, see © and ©, while (C) respectively 
gn}m (C) are defined, in terms of the quantities fn ^ (C) respectively (C) (see © respectively ©), as follows: 


f(k) (C) = 9 f ( n ] (Z) 

Jn,m \S_J o 


Jj) __ d 9 { n fe) 


qd' (c\ = 
1 yn,m V_1 / 


£=C 


2=C 


(9) 


Expressions of fn)m (C) and On In (C) w ith j = 1, 2,3,... are reported in the Appendix, see (l64l) and (1651) . 

Then the N eigenvalues A m of the N x N matrix L (C), 

L (C) (C) = A m v ("*) (C) , m= 1,..., N (10a) 

hence the N roots A m of the following polynomial equation (of degree N in A): 

det [L (C) - A] = 0 (10b) 

—are given by the formula 

A m (/3i, P q ) = m fl (/3fe — 1 + to) , m = l,...,N .D (10c) 


Remark 2.2. The functions fn.ln (C) and gn]m (C) are universal: they do not depend on the generalized 
hypergeometric polynomial under consideration (see Remark 2.1). But of course their arguments do, being 
the N zeros ( n of the polynomial P N (aq,..., a p \ /3 1 , ..., /3 ? ; £) , see ©. □ 

Remark 2.3. The N eigenvalues \ m of the N x N matrix L (see ©) depend only on the q parameters 
P k (see (nodi '), while the matrix L depends itself on the q + p parameters /3 fe and aj —via the dependence of 
the parameters b k respectively aj on f3 k respectively aj (see © respectively ©) and the dependence of the 
N zeros C, n on the parameters /3 k and aj, see (12bl) or, equivalently, on b k and aj, see ©. Hence the N x N 
matrix L is isospectral for variations of the p parameters aj. And note moreover that the N eigenvalues A m are 
integer (or rational) numbers if the q parameters j3 k are themselves integer (or rational) numbers: a nontrivial 
Diophantine property of the TV x N matrix L. □ 

Remark 2-4- All the above results are of course true as written only provided the N zeros ( n are all different 
among themselves; but they clearly remain valid by taking appropriate limits whenever this restriction does 
not hold. □ 

Remark 2.5. Immediate generalizations—whose explicit formulations can be left to the interested reader— 
of Propositions 2.1 and 2.2 obtain from these two propositions via the special assignment aq+j = /3p+j 
for j = 1with r an arbitrary nonnegative integer such that both q = q — r and p = p — r are pos¬ 
itive integers. These propositions refer then to the N zeros of the generalized hypergeometric polynomial 
Pjv (oq, ap\ ..., (3^, z )—which depend only on the 1 + p + q = 1 + p + q — 2r parameters N, aj with 
j = 1, ...,p = p — r and b k with k = 1 , ...,q = q — r, but feature quantities b k and aj (see © and ©) that 
depend on the 1 + p + q parameters N, aj with j = 1, ...,p and b k with k = 1,..., q. □ 

The two Propositions 2.1 and 2.2 are proven in the following Section 3; some comments and prospects of 
future developments are outlined in the last Section 4. 

Let us end this Section 2 by displaying explicitly the above results for small values of the integers p, q and 
(of course) r (see Remark 2.5). 


2.1 The case p = q = 1, r = 0 

For p = q = l, r — 0 (for the definition of r see Remark 2.5) implying (see © and ©) 


a 0 = aq, a 1 = -1, b± = /3 X - 1, b 2 = -1 , 


let the N numbers £ n be the N zeros of the hypergeometric polynomial 

N 

P N (aq -,Pp,z) = ^2 

m =0 


(fil )m 


( 11 ) 


( 12 ) 
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Then Proposition 2.1 implies that these N zeros C >n satisfy the following system of N nonlinear algebraic 
equations 

N-l-a 1 +P 1 C„-2 (C„-l) <4 M) (C)=0, n = 1, N . (13) 

Notation 2.2. Above and hereafter the notation a n’^ (C) is defined by (1521) . □ 

As for Proposition 2.2, it implies in this case that the following N x N matrix L ((f) , defined componentwise 
as follows (see © with (1M1) . (RT51) and (I52H . 


N 


Lnm (C) — &nm i /3i + 2 £ 


Ce (Ce — 1) 


-2 (1 - 6 nm ) 


(Cn - Ct) 
Cn (C„ - 1) 


features the N eigenvalues 


(C -C ) 

VSn Sm/ 

A m = m (Pi — l+m) , m=l,...,N. 


(14a) 

(14b) 


2.2 The case p = 2, q = l, r = 0 

For p = 2,q = l,r = 0 (for the definition of r see Remark 2.5), implying (see © and ©) 

a 0 = op a 2 , ai = - (op + a 2 ), a 2 = 1 , 6i = /3 X - 1, b 2 = -1, (15) 


let the iV numbers be the iV zeros of the hypergeometric polynomial 


N 

Pn {a 1 ,a 2 \P 1 ;z) = 

m =0 


' RU a l)m ( Q 2) m ^ 

(^l)m 


(16) 


Then Proposition 2.1 implies that these AT zeros C tn satisfy the following system of N nonlinear algebraic 
equations 


- oq a 2 + {N - 1) (qi +a 2 + 1) +/3 1 ( n 
+2 (3 - N + a\ + a 2 — C„) Cn' 1 ^ (C) + 3 &n’ 2 \() + 3 (C)] 2 = 0 , 

n = 1, ..., N . (17) 

As for Proposition 2.2, it implies in this case that the following N x N matrix L (£), defined componentwise 
as follows (see © with HMD, (ESI) and ([52]) ). 


Lnn (C) 


& + (JV - 3 - a, - «,) „<y> (<) + <ri=.=> (<) - 3 <r<«> (f) 


+3 4 1 ’ 1 ' (<) s','- 21 (<) , n = 1,2,.... N 


(18a) 


(C) = 2C„ 


+ «2 — N — £ n — 3 ah ’ "* (C) 


3 Cn 


(C — C r 

\Sn Sra/ 


(Cn - Cm) 

, n,m=l,...,N, n^m 


features (again) the N eigenvalues 


A m = m (/3 1 — l + m) , m = l,...,N . 


(18b) 


(18c) 


Note the isospectral character of this matrix L (C, ) , which depends explicitly on the 2 parameters a = ai+a 2 
and /3 1 and implicitly on the 3 parameters a\,a 2 and P 1 via the dependence on these 3 parameters of the N 
zeros ( n of the polynomial Pn (op, a 2 ] /3p,z), while its eigenvalues A m only depend on the single parameter /3 1 . 
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2.3 The case p = q = 2, r = 0 

For p = g = 2,r = 0 (for the definition of r see Remark 2.5), implying (see (0 and 0) 


a o — ot\ a 2 , cli — — (ai + a 2 ) , 02 — 1 

C = (1 — /Si) (1 — 0 2 ) , = 2 /5 2 


= 1 


let the N numbers ( n be the N zeros of the hypergeometric polynomial 

N 

P N (op ,a 2 ;P 1> P 2 \z) = Y 


N r (-N) m ( ai ) m (a 2 ) m /"“I 


m —0 L 


m ! (/3i) m (/3 2 ) r 


(19) 


( 20 ) 


Then Proposition 2.1 implies that these N zeros C, n satisfy the following system of N nonlinear algebraic 
equations 


— ai« 2 + (N — 1) (an + a 2 + 1) + 0i 0 2 Cn 
— 2 [(1 + 01 + /3 2 )C n — OL\ — OL2 + N — 3] cr^ 1 ’ 1 ^) 


+3 (C„ - 1) i cr^CC) ~ a n ' 2> (() r' = 0 , n=l,...,7V. 


1 2 


( 2 , 2 )/ 


( 21 ) 


As for Proposition 2.2, it implies in this case that the following N x N matrix L (C), defined componentwise 
as follows (see 0 with (l64l) . (1651) and (1521) 1. 


+2 


Lnn (c) = 01 02+ [ 5 + 2 (0! + 0 2 )} ^ (C) + 6 a™ (C) 
(N — 3 — ai — a 2 ) a£’ 2) (C) - 3 o£’ 3) (C) + 3 oi M) (C) c#- 2) (C) 


-3^ 1} (C) k 1 - 1 ) (C) + 2 (C) 


n = 1,2,...,7V 


( 22 a) 


(C) = 2 Cn 


ai + a 2 — TV + (2 — /3 1 — /3 2 ) Cn + 3 (Cn — 1 ) ah ' (C) 


3 Cn (Cn ~ 1) 
(Cn - Cn) 3 


(Cn - Cm) 

, n,m=l,...,N, n^m 


features the N eigenvalues 


A m = m (/^ — 1 + m) (f3 2 — l + m) , m=l,...,N. 


( 22 b) 


( 22 c) 


Note the isospectral character of this matrix L (C) , which depends explicitly on the 3 parameters a = 
ai + a 2 , 0 1 and /? 2 and implicitly on the 4 parameters ai,a 2 , 0 1 and /? 2 via the dependence on these 4 
parameters of the N zeros C n of the polynomial Pn (oq, o 2 ; 0 1 ,0 2 , z), while its eigenvalues A m only depend on 
the 2 parameter /3 X and /3 2 . 


2.4 The case p = q = 2, r = 1 

For p = q = 2, r = 1 (for the definition of r see Remark 2.5) we have (as above, see (TT9l) . but now with 
02 = “ 2 ) 


a o — ai a 2 , ai — — (ai + a 2 ), a 2 — 1 , 

bi = (1—0^) (1 — a 2 ) , b 2 =2 — /3 1 —cx2, 63 = 1 , (23) 

and we now define the N numbers Cn as the N zeros of the polynomial (fl2l) (rather than (121)1) : so that these 
N zeros do not depend on the arbitrary parameter j3 2 = a 2 ). Then Proposition 2.1 implies that these N zeros 
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C n satisfy the following system of N nonlinear algebraic equations 


- ai a 2 + (N - 1) (ct\ + a 2 + 1) + /3i a 2 Cn 
+2 [3 — N + a.\ + a 2 — (1 + Pi + a 2 ) Cn] ^n’^lC) 


+ 3(Cn - 1) 


r(M)/ 


(0 - \ = 0 , n= 1 ,..., N 


(24a) 


—which is different from m, although satisfied by the same zeros. But since this system of N equations 
must hold for arbitrary values of the parameters a 2 , it amounts in fact to the following two separate systems 
of N equations: 

- ai + N- 1 + /3, Cn + 2 (1-CJ ^ M) (O = 0 , n = 1 ,N , (25a) 


(N- 1) (o! + 1)+2 [3-N + a 1 -(l + p 1 ) C„] ^i M) (0 


+ 3(Cn _ 1) 


T"( 151 ) t 


(C) -^ 2 ’ 2) ( C))=0, n=l,...,N 


(25b) 


which must both be satisfied by the N zeros C n of the polynomial m- Indeed the first of these two systems 
coincides with (fl3l) ; while the second is new (but in fact both these systems of N equations are not quite new, 
see Subsection 2.5). 

As for Proposition 2.2, it implies in this case that the matrix L (£) defined by (I22al) and (122bl) (of course 
with /J 2 = a 2 ) hence reading 


+2 


Lnn (C) — Pi 012 + [5 + 2 (/3 1 + a 2 )\ O’n’ 2 ' 1 (c) + 6 cr^ 3,3) (c) 
(N - 3 - ar - a 2 ) a^’ 2 ) (C) - 3 a^’ 3 ) (C) + 3 a^’ 1 ) @ @ 


-3.1 U) (C) k 1 - 1 ) (C) + 2 a( 2 > 2 > (C) 


n = 1,2,..., AT 


(26a) 


(C) = 2C„ 


ai + a 2 — N + (2 — P 1 — a 2 ) ( n + 3 (C„ — 1) cri 1 ’^ (C) 


3 Cn (Cn - 1) 


(Cn - c m r 


(Cn - Cm) 

n, m = 1,..., N , n m 


(26b) 


—but now with the N zeros ( n in these two formulas being again those of the generalized hypergeometric 
polynomial m rather than (BUI) , so that they only depend on N, ai and (but not on /? 2 = a 2 )—features 
the N eigenvalues 

A m = m (/3 1 — l + m) (a 2 — 1 + m) , m=l,...,N. (26c) 

But again, since these properties must hold for arbitrary values of the parameter a 2 , they amount to two 
separate statements, the first of which is easily seen to reproduce the statement that the matrix (lldal) features 
the eigenvalues (I14bl) . while the second states that the above matrix (I26al) . (I26bl) with a 2 = 0, i. e. 


+2 


Lnn (C) = (5 + 2 ft) ai 2 ’ 2 ) (C) + 6 ai 3 ’ 3 ) (C) 
(N- 3 - ai ) ai 1 ’ 2 ) (C) ^ 3 ai 2 ’ 3 ) (C) + 3 ai 1 ’ 1 ) (C) (C) 


-3 4 14) (C) k’ 1} (c) + 2 ai 2 ’ 2) (C) 


n = l,2,...,N 


(27a) 


Lnm (C) — 2 C„ 


ai — N + (2 — ft) Cn + 3 (Cn - 1) ^ (C) 


3 Cn (Cn - 1) 


(Cn - Cm) 


(Cn - Cm) 

n,m=l,...,N, n^m 


(27b) 
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features the eigenvalues (126el) with a 2 = 0 , i. e. 


A m = to (to — 1) (/3 1 — 1 + to) , m = l,...,N. (27c) 

Note that, while these eigenvalues depend on the parameters /3 1 , they do not depend on the parameter 
op; hence in this case the matrix L (£), which itself depends on the 2 parameters op and /3 1 , is isospectral for 
variations of the parameter a\. 


2.5 Results for Jacobi polynomials 

For p = q = 1—or, equivalently, for p = q = 1 (for this notation see Remark 2.5 )—the generalized hypergeo¬ 
metric polynomial is simply related to the Jacobi polynomial (see eq. 10.8(16) of [16) J: the transformation (up 
to an irrelevant multiplicative constant) from the generalized hypergeometric polynomial Pat (cm; /3 X ; z), see 
m, to the standard Jacobi polynomial Pjf’^ ( x) (see [16]) corresponds to the change of variables /3 X = a + 1, 
ol\ = N + a + /3 + 1 and z = 2/ (l — x). It can thereby be verified (with some labor) that the results— 
corresponding to Proposition 2.1 with p = q = 1, r = 0 respectively p = q = 2, r = 1—reported above for 
these cases reproduce known results [2J: specifically (fUl) (or, equivalently, (12 5 all ) respectively (I25bl) reproduce 
(up to appropriate notational changes) eqs. (5.2a) respectively (5.2b) of [2]. 

On the other hand Proposition 2.2 with the above assignments of p , q and r seem to produce new results 
for the N zeros x n of the Jacobi polynomial Pjf’^ ( x ), as displayed below. Indeed for p = q = 1, r = 0 it 
yields the following 

Proposition 2.5.1. The N x N matrix L ( x) defined componentwise, in terms of the N zeros x n of the Jacobi 
polynomial P^*’^ ( x ) and the two parameters a and (3 as follows, 


N 


Lnm (ifi) — $nm ^ O J 1 J ^ ] 

1=1, i^n 


( 1 + 2 ~t) (1 — x n ) 


(x n - Xc) 


- (1 - Sum) 


(1 + X n ) (1 X m ) 

(Xn Xrn) 


(28a) 


has the N eigenvalues 

A m = to (to + a) , to = 1 ,..., N . □ (28b) 

Let us again note the isospectral property of this matrix L(x ), whose elements depend, via the N zeros 
x n = x n (a, /?), on the two parameters a and /3, while its eigenvalues depend only on the parameter a. 

Also note that Corollary 5.2.2 of j2] with s = N — l,r = 0, and of course n = N, identifies an N x N matrix 
G (see eq. (5.17) of [5] with C and X defined by eq. (5.4) and (5.3) of [2]) defined componentwise as follows: 


G n m — 5 


N 

nm ^ ( 
1=1 


(1 - x}){l ~ X n ) 

(X n - Xt) 2 


- (1 - $nm) 


(1 Xrn) (1 + Xm) 
(Xn Xm) 2 " 


(29a) 


and states that its N eigenvalues g m read as follows (see eq. (5.19) of [2], with N — to replaced by to — 1 since 
these numbers span the same set of values—from 0 to N — 1 —for to = 1,..., N): 


gm = (to 1) (to + a — 1) , TO = 1,..., TV . 


(29b) 


These formulas, (l29l) . are similar to, but different from, (|28[> : although of course there must be a way to relate 
them, since they hold for the same set of N + 1 numbers a and x n = x n [a, /?). 

And likewise for p = q = 2, r = 1 it yields- via a treatment analogous to that of Section 2.4—the following 
Proposition 2.5.2. The N x N matrix L (x) defined componentwise, in terms of the N zeros x n of the Jacobi 
polynomial P^'^ (x) and the two parameters a and (3 as follows, 


1 ■ r 


= (7 + 2a) c>A 2 ’ 2) (x) + 6 cr ( n 3 ’ 3) (x) - 2 (4 + a + (3) a ^ (x) - 6 cr^ 2 ’ 3 ) (x) 


1 2 


<r { n’ 1) (x) - 6 a£’ 1 ) (x) (®) + 6 cr^’i) 0 - 6 - 2 ) ^ ? 


-3 


(30a) 
















(a + /3 + 1) (1 - x n ) + 2 (1 - a) + 3 (1 + x n ) cr^ 1 ’^ (x) 


1 - Xr, 


L, n •A'm 


-3 (l + x n ) 


1 %m 
%n %m 


where (see ([52]) ) 


N 


r (r, P ) {3 . )== £ 


1 - Xg 


r—p 




has the N eigenvalues 


£=1 -i^n 

X m = m (to — 1) (to + a) , to = 1 ,N . □ 


(30b) 

(30c) 

(30d) 


Note again the isospectral character of this matrix, which depends on the two parameters a and (3, while its 
eigenvalues depend only on the parameter a. 

3 Proofs 

In this section we prove the findings reported in the preceding Section 2. 

Let the t-dependent monic polynomial, of degree N in z and characterized by its N zeros z n (t) and its N 
coefficients c m (t ), 


JV 


N 


VhvOM) = n [z-z n (t)} = z N + Z m ] 

" =1 rrt—1 

N 

= ^2 [Cm ( t ) z N ~ m ] with Co = 1 , 
m =0 

satisfy the linear Partial Differential Equation (PDE) 


(31) 




dt 


d \ p 

d,n 


d z 


aj - I z - - N 

o z 


3 =1 

Tpl vOM) ■ 




(32) 


Clearly this implies that its coefficients c m ( t ) satisfy the following system of N linear Ordinary Differential 
Equations (ODEs): 

d N 

d t . 


) C m (t) = TO 

n (Pj - 1 + m) 

/ 

i=1 


+ (N + 1 — to) 

with the conditions (see dTTTl) 1 


n ( a j -1+to) 

j =i 


Cm (t) 

Cm-1 (t) , TO = 1, ..., N , 


C 0 (t) = 1 


c j {t) = 0 for j < 0 and for j > N . 
Hence the general solution of this system of N linear ODEs reads as follows: 


N 

c(f)= £ [dm exp (X m tj 

m= 1 


-(m) 


(33a) 

(33b) 

(33c) 

(34) 
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where the N (t-independent) parameters rj m can be arbitrarily assigned (or adjusted to satisfy the N ini¬ 
tial conditions c m (0)), while the numbers A m respectively the N -vectors v^ m ' 1 are clearly the N eigenvalues 
respectively the N eigenvectors of the algebraic eigenvalue problem 

A u( m) = A m fj( m ) , m = l,...,N, (35a) 

with the N x N matrix A defined componentwise as follows: 

A m,m = m Yl (Pj -1 + m) , m = 1 ,N , (35b) 

i=i 


P 

A m , m -1 = {N + 1 - to) n (aj-l + m) , m = 2,...,N, (35c) 

j =i 

with all other elements vanishing, Am,n = 0 unless n=morn = m-l. 

The (lower) triangular character of the matrix A implies that its N eigenvalues X m can be explicitly 
evaluated: 

<7 

Am = A m,m=m T\ (Pj -1 + ITl) , TO = 1, N . (36) 

3 = 1 

Let us now denote as 4> N (2) the t- independent monic polynomial solution of (15^1) . 


N JL 

i’N (z) = U[z- U =**+£ [7m Z N ~ m ] 


n =1 
N 


m =1 


= bm Z N m ] with 70 = 1 ) 


(37) 


771=0 


hence the monic polynomial solution of the ODE 


‘-:- N ) jJi L 




-(A) h 


d 2 


j=i L 




1p N {z) = 0 • 


(38) 


Note that we denote as its coefficients and as C >n its N zeros, see (l37l) (the fact that the notation for the zeros 
is identical to that used above, see (I2b[) . is not accidental: see below). Clearly the t -independent coefficients 
7 m correspond to the ’’equilibrium configuration” c m ( t ) = of the linear ’’dynamical system” (133al) . hence 
they are characterized as the solutions of the system of N linear algebraic equations 



9 . 


V 

(to + 1) 

FI (/3j + m) 

3= 1 

7m-}-1 = (TO - N ) 

n ( aj +to) 

3= 1 


(39a) 


with the conditions 

7o = 1 . 7iv+i = 0 • (39b) 

The first of these conditions corresponds to (I33bl) : the second (which via the recursion II 39 all clearly implies 
7 m = 0 for m > N) corresponds to (133cl) and it is automatically satisfied because the right-hand side of the 
recursion (139a|l vanishes for to = N due to the factor (to — N) (and note that no conditions need to be assigned 
on 7 m with to < 0 since no such values enter in the recursion (I39al) ; in any case for to = — 1 the recursion (139al) 
would imply ')_ 1 = 0 since its left-hand side vanishes due to the factor (to + 1 )). 

It is then plain from the two-term (hence explicitly solvable) recursion relation (l39l) that the 7V+1 parameters 
7 m read as follows: 

( N )m n (aij) m 

7m =- ^ -> rn = 0,1,..., N. (40) 

m! n {Pi) m 

3 =1 
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This—besides implying that the ODE (1551) does possess a polynomial solution ip ( z ) of degree N in 2 —shows 
that this f-independent polynomial solution ip (z) of the PDE (1321) coincides with the generalized hypergeometric 
polynomial (l2all (compare (f2aT) to (1571) with ([40]) ); of course up to an overall multiplicative constant, which 
can be arbitrarily assigned due to the linear character of the ODE (1551) . and was chosen above so that the 
polynomial ip (z) be monic, see (EH, hence indeed coincide with the generalized hypergeometric polynomial 


Our next task is to identify the ’’equations of motion” characterizing the t-evolution of the N zeros z n (t ) 
of ip (z , t) , see (f3Tlt . implied by the PDE (1551) and by the corresponding f-evolution of the N coefficients c m ( t ), 
see (1331) . To this end it is convenient to firstly reformulate the PDE (l32|) as follows: 


(§i) ’ M + ) 



Vhv(M) , (41) 


where of course the new parameters bk respectively aj are related to the parameters /3 k respectively otj as 
detailed above, see and ©. 

Then it is easily seen that the equations of motion characterizing the f-evolution of the N zeros z n (t) of 
ip ( z,t ), see (15T1) . read as follows (of course below a superimposed dot denotes a ^-differentiation): 


9+1 


Zn = 


fc= 1 


Y [ bk fn ] U) - Y [ a J 9n ] U) 

3 =0 


n = 1 ,N 


(42) 


Indeed, this clearly follows from the PDE (TUI) via the identities (I58al) and (160al) and the analogous identity 
corresponding to the logarithmic f-derivative of (13TT) hence reading (via the short-hand notation (1151) ) 



V’jv ( z ,t) •$=> ~Zn ■ 


(43) 


It is moreover plain from the developments reported above (see (1571) and the sentence following (l40l) ) that 
the coordinates C, n characterizing the equilibrium configuration £ = (Ci) •••) Cn) °f this system, —which of course 
(see (1151) 1 satisfy the set of N algebraic (generally nonlinear ) equations ©— coincides with the N zeros of the 
generalized hypergeometric polynomial Pn (oq,..., a p ; j3 1; ..., /3 q \ z), see (I2bl) . 

Proposition 2.1 is thereby proven. 

Our next step is to consider the behavior of the dynamical system (1421) in the infinitesimal vicinity of its 
equilibrium configuration z ( t ) = £. To this end we set 

z{t)=( + ex(t); z n (t) = C„ + £ x n (t) , n = l,..,N , (44) 

with e infinitesimal. We thereby linearize the equations of motion (1421) . getting 

N 

X_ — I X , X n — ''p ) ( Lnm *£m) 5 Cl 1, ..., IV , (45) 

771=1 


with the N x N matrix L defined componentwise by ©. 

The general solution of the system of linear ODEs (ll5l) reads of course then as follows: 


N 

x(t)=Y [Vm ex P t) H (m) 

m =1 


(46) 


where the N (t- independent) parameters can be arbitrarily assigned (or adjusted to satisfy the N initial 
conditions x n (0)), while the numbers X m respectively the IV-vectors are clearly the N eigenvalues respec¬ 
tively the N eigenvectors of the algebraic eigenvalue problem (llOal) . But the behavior of the dynamical system 
m in the immediate vicinity of its equilibria cannot differ from its general behavior, which is characterized 

by the N exponentials exp ^A m t'j , as implied by the relation between the N zeros z n (t) and the coefficients 
c m ( t ) of the monic polynomial (of degree N in z) ip N ( z,t ), see (1511) . and by the explicit formula (1M1) with 
(1551) detailing the time evolution of the N coefficients c m ( t ). Hence the (set of) eigenvalues A m of the matrix 
L, see ©, must coincide with the (set of) eigenvalues A m , see (l36l) . of the matrix A, see (l35l) . 

Proposition 2.2 is thereby proven. 
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4 Outlook 


A follow-up to the present paper shall take advantage of the new polynomial identities reported in Appendix 
A in order to identify new classes of solvable IV-body problems of ” goldfish” type mmmm involving 
many-body interactions and featuring several free parameters, by extending in a fairly obvious manner (indeed, 
see (01) the approach and findings reported in US- 

Two possible directions of further investigation shall try and extend the approach and findings, reported in 
this paper for the N zeros of hypergeometric polynomials of order N, to the N zeros of basic hypergeometric 
polynomials of order N and to the, generally infinite, zeros of (nonpolynomial) generalized hypergeonretric 
functions. 
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6 Appendix A: Polynomial identities 

Several polynomial identities are reported in Appendix A of |14j (see in particular the paperback version of this 
monograph, where the identities (A.8k) and (A.81) are corrected) and in Appendix A of [15]. In this appendix 
we report, and then prove, two additional classes of analogous identities. These findings have an interest of 
their own: for instance they open the way to the identification of new classes of solvable IV-body problems of 
’’goldfish” type (as mentioned in Section 4). 

To make this Appendix self-consistent we firstly introduce the notation we use below, even if this might 
entail a bit of repetition (see [14]). 

Let if) ( z ) be a monic polynomial of degree N in z, and denote as z n its N zeros, 


N 

( z ) = II ( z ~ Z n ) , 

n =1 

and as z the N -vector of components z n , z = (zi,Z n). 

We then introduce the following short-hand notation: the formula 


stands for the identity 


D 


D V> {z) f n (z) 


N 

v* (*) = £ [( z - z «) _1 a te) 


71=1 


Here D is a differential operator with polynomial coefficients acting on the variable z, say 

J r 


3 = 0 


(z) 


d 

d z 


with the coefficients (z) polynomial in z. The simpler examples of such identities are 


(47) 

(48a) 

(48b) 

(48c) 

(49a) 
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(49b) 


Z d~z~ N ) 

clearly the first obtains by logarithmic differentiation of (H71) . and the second from the first via the identity 
z/ (z — z„) = 1 + z„ / (z — z„) ("these identities. (I49al) andl49bll. are reDorted in as (A.4) and (A.6a)). 

Finally, to write more compactly some of the formulas obtained below and reported in Section 2 we introduce 
the following convenient notations: 

N N 

= V 


E - 




1=1 ;(n) £=l\£±n 

N N 

E - E 

£,m=l\(n) 

N N 

E - E 

£,k,m=l’,(n) £,k,m=l;£^n,k^n,m^n,£-^:k,£^m,k^m 


(o\ 

d n{tm) U) = ( z n ~ Zg) ( z n ~ Z m ) (Zg - Z m ) 


d (i) 

U 'n(^fcm) 


(z) = {z n - Zt) ( Zn - Z k ) (z n - Z m ) {zg ~ Z m ) (zg - Z k ) ( Z k - Z m ) 


(50a) 

(50b) 

(50c) 

(51a) 

(51b) 


Remark A.l. The upper label on the quantities d^ m ) G=0 an d ^n(gkm) (-) a reminder that (l) 

(z) involves the coordinates of the 4 


<4) 


j(3) 


involves the coordinates of the 3 different zeros z n . zg, z m , and dt 


(4) 


i(£kr 


different zeros z n , zg, z k , z m . It is moreover plain that these quantities—that shall play the role of denominators 
in sums over the zeros, see below—are antisymmetric under the exchange of any pair of these zeros: this entails 
that, by taking advantage of these (anti)symmetries, the arguments of the sums reported below could be 
rewritten in many different ways by appropriate exchanges of the dummy indices being summed upon. We use 
this property to also obtain convenient expressions based on the quantities (Jn' p ^ (z), introduced immediately 
below. □ 

It is now convenient to introduce the quantities 


r ( n r ’ P) U) = 


N 

E 

g=l-g^n 


_(z n - Zg) P _ 


(52) 


Here and hereafter r and p are arbitrary nonnegative integers (unless otherwise indicated). Clearly this defini¬ 
tion implies 


d a 


C.p) 


(!) 


= - oS < 7 (r, P+ l), ) (1 _ § s [r Zn + (p- r) Zm] Z* 

— P Onm O n t!J T Onm) \n+l 

[.Zn ^ Z m ) 


hence 


d On’ r) (z) 
d z m 


= -r 6 nm <J^ r+1) (z) + (1 - 6 nm ) 


(Zn Zm) 


r+1 


and it is also obvious (and used below) that 


^°’ 0) (z) = N- 1 . 


(53a) 


(53b) 


(54) 


Remark A.2. The quantities crf’^ (z) , see HS2D, are related to each other by trivial identities, for instance 
clearly the replacement of z\ with zl -1 [z n + (zg — z n )] in the numerator of their definition (1521) implies 


^ { n ,p) (z) = z n hp) (z) - <r£ 1>p 15 (z) , r,p = 1,2,3,... 


(55) 


yielding (by iteration) 


rt p) (!) = E 

k—0 


(-if 


s-fc (r-s,p-k) 


(!) 


s = 0,1,2,..., min (r,p) , 


(56a) 
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hence in particular 


r (r,p) ( 7 \ = 


fe) = E 


,(r,p) \ _ 


k =0 
P V 


<-+ (I) <r " ■ U) 


fe) = E 


k—0 


(~l) k ( P \ z p ~ k ak-P’P-V 


GO 


if r < p 


if r > p . □ 


To obtain the second version of the identities reported below we also used the identity 

ri’ p) GO = <r£ +1 ’ri (z) + a^’P-V (z) , r = 0 , 1 , 2 ,.... p = 1 , 2 ,3,... 


(which correspond clearly to (l55l) h and the property ( 1 + 0 ) . 

The first class of new identities reads as follows: 

^ ( Z ) fn 3) fe) , J = 1) 2 , ... , 

with the TV quantities fn ■* (z) defined recursively as follows: 


with 

implying 


N 


fn +1) GO = - fn ] U) + E 

t=l;(n) 


Zn (z) + Z t f^ J) (Z) 


(J) 


Z n ~ Zl 


fn ] U) = Z n 


N 


= 


fn ] GO = -Zn + 2 E 

£—1 ;(n) 

-1 + 2 (z) 


Z n Zj 
Zn - Zt 


3 = 1 , 2 ,... 


AT 


/« 3) GO = Z„ - 6 E 


Z n Zt 


N 


+ 6 E 


Z n Zg Zqjl 


^ Zn ~ Zp , 

£—1 ;(n) v 7 ^,m=l;(n) \ u, n(^' 


EE GO 


= < 1 - 6 cr£’ 1) (z) - 3 c>i 2 ’ 2) (z) + 3 lo 41 ’ 13 CO 


1 2 


AT 


/n 4) GO = + 14 E 

£=1 ;(n) 
N 

+24 E 


Zl 

Zn - ^ 


-36 E 


N / o ' 
Zn z| Z m 


*■ —' \ j( 3 ) ( r \ 

l,m=l;(n) \ u ntm \^> 


Zn Zt Zl Z ^ 


£,k,m=l-,(n ) \ u n£km V—/ / 

= |-1 + 14 (z) + 18 cri 2,2) U) + 8 <E’ 3) (0 

GO] - 12 E’ 1} (0 E’ 2) GO + 4 \ (J n ,1) GO 


-18 


and so on (via standard, but increasingly tedious, computations). 
The second class of new identities reads as follows: 



31 

EEO 

1 


, , i 2 -J- ~ n) 00 9n ] GO , j = 1, 2,... 

a z J \ a z ) 


(56b) 

(56c) 

(57) 

(58a) 

(58b) 

(59a) 

(59b) 

(59c) 


(59d) 


(60a) 
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(for the case j = 0 see (I49al) implying g(z) = 1), with the N quantities (z) defined as follows: 


N 


9; 


(i) u)= E 


ti j) go+ n j> go 


O') 


Z n - Zt 


t=l ;(n) 

implying (for the second version of these formulas via ([54]) ) 

N 

9n ] U) = Zn E 

i=l ;(n) 

= iV-l+2 f 7£ 1 - 1 )U) , 


, j = 1,2,... , 


1 


Z n ~ Zt 


N 

E 

t=l;(n) 


Zt 


Zn ~ Zt 


(60b) 


(61a) 


N 


,( 2 ) 


£—1 ;(n) 


( S )- 2 E 


AT 


2 m + 


u) 


= 1 - TV + 2 (N - 3) a*, 1 - 1 ’ (z) - 3 a^ 2 ’ 2 ) (z) + 3 GO 


AT 


£=1 ;(n) 


>(?>fe> = £ +« y 


- 2 / 


N 


£,m=l\{n) 


(Z n + Zf) zl 

</ 3 A > ( 0 ) 

n(£m) v —/ 


AT 


(Zn + Zl) 2 ^ z% 

d nltkm) U) 


+6 E 

£,k,m=l](n) 

= N - 1 - 2 (3 N - 7) a™ (z) - J (IV - 7) ^ 2 ’ 2 > (z) + 6 (z) 


+ 4 (JV-7) 


+ s Ikk) 


(61b) 


(61c) 


and so on (via standard, but increasingly tedious, computations). 

To prove the identities (l58l) we note first of all that (15 9 all coincides with (149bl) (via (|48l) h Next—to prove 
(|58al) with (I58bl) —we apply the operator (z d/dz — N) to (the long-hand version of) (|58al) . getting thereby the 
following chain of equations: 


j+1 


( z - n) [${z) E ~ Zn ) 1 f™ 3) fe) 

' ' f n—1 

N 

E [(z - Z„) _1 U) 


(62a) 


Z __AT] ^(z) 


J n=l 


-V 7 (2) E [ z ( 2 - ~«) 2 k fe) 

n=l 

^ {z) | E [(- ~ ~ n ) _1 _ ~<o _1 ^ & 

^ n,£=l 

AT ^ 

“E [~ ("-^)~ 2 fn j) GO 


(62b) 


n=l 


(62c) 
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V’ (z) | ( Z ~ z ™> 1 ( 2 “ 1 ^ fn } fe) 

- 51 [(* _ Z ") _1 fn ] fe) 

n=l 


(62d) 


V> (z) j 5Z (* “ ~«) 1 “ z i) 1 ^ fn j) U) 

^ti,&= 1 ; £ 7^71 
AT 

~ 5 Z (* “ ^) _1 _ ^) _1 ^ fn ] U ) 

n,£=l; £^n 


N 1 

[( 2 - 2 «) _1 fi 3) ( z ) > 

n=l J 


(62e) 


f ^ 

V >( z ) s 5 Z [( 2_2 ") _1 {'/ n ^ U ) 


<71=1 


N 


+ 5Z (*« - z t) 1 

£=1; 


^ fn ] U) + Z n fl 3 ’ (z) 


U) 


(62f) 


The first equality is implied by (the long-hand version of) (|58ajl : the second by standard differentiation; 
the third, by (the long-hand version of) (I49bl) : the fourth obtains by taking advantage of the cancellation 
of the terms featuring double poles; the fifth, by using the elementary identity (z — z n ) 1 (z — zi) 1 = 
{z — z n ) 1 — (z — ze) 1 [z n — zt) ; and the sixth and last, by exchanging the two dummy indices n and i 

in the second sum. The last equality is clearly the long-hand version of (158aD (with j replaced by j + 1 and 
with (I58bl) ): the recursion (I58bl) is thereby proven. 

To prove (IfiOap with (IfiObl) we apply the operator d/dz to (the long-hand version of) (|58ajl . getting thereby 
the following chain of equations: 


_d_ 

d z 

_d_ 

d z 


d V 

2 l~z ~ N ) ^ ^ 
N 


,i>( Z ) 51 “ Z ") 1 fn } U) 


(63a) 


ip(z) 


N 

5] \(Z - ZnY 1 fY (Z) 


J 71=1 


N 

-V’ i z ) 51 [(~ _ Z n )~ 2 fn j) U) 


71=1 


(63b) 


i’(z) { ( z ~ z ^ 1 ( z ~ z t ) 1 fnHz)} 

71, £=1 

N 

( Z ) 51 [( Z ~ Z ^~ 2 fn' 1 (z) 


71=1 
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(63c) 






































N 



= i>{z) 

n 

E (^ “ Zn ) 1 ( Z “ Z ^> 1 fn ' U) 

4=1; l^n 

(63d) 

V’(-) 

N 

E 

n,^=l; £^n 

| [(z - Z„) 1 - (^ - Zl) 1 

( Zn-Zt ) 1 /1 J) U)} 

(63e) 

■0(2) 

N 

E 

|( 2 - 2 :„) _1 (z„ - zi)~ l 

U) + (*)] } • 

(63f) 


n,£=l; £^n 


This sequence, HMD, of equalities is sufficiently analogous to that reported above, see (15^1) . not to require any 
repetition of the explanations provided above (after (1C>21) 1. And clearly the last equality proves the identity 
(160al) with (I60bl) . 

Finally, let us report explicit expressions of the quantities fn,m (C) and 9n}m (C) , see © (conveniently 
obtained via (l53l) from the expressions of fn' 1 (£) and (C) i n terms of the quantities dn’^ (C), see (l59l) and 

HD): 


fk% (c) = 5 


nm ■> 


(64a) 



Sr, 


1 + 2 a. 


( 2 , 2 ) 


(c) 


+ (1 — S nm ) 2 



(64b) 


A%(0= 4™ {l + 9^»(S)+I^(£) 

-3 ■» (C) [+ 1 ) (<) + 2 (<)] } 


-6 (1 - Snm) 



(Cn - Cm) ^’ 1) (C) 
(Cn - Cm)' 


and so on (via standard, but increasingly tedious, computations); respectively (see © and (l60l) l 


9 ( nl (C) = -2 5 nm (C) + (1 - S nm ) 


2 Cn 

(Cn - Cm)' ’ 


(64c) 


(65a) 


(C) = Snm {-2 (N - 3) (C) + 6 (C) 

-6 ^ (C) ^’ 2) (C)} 


+ (1 


$nm) 


2 Cn [(N - 3) Cn - N Cm] 

(Cn - Cm) 3 


6 Cn 0'»’ 1) (C) \ 
(Cn Cm) 2 J ’ 


(65b) 
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sgl (C) = inn, {2 (3 N - 7) (C) 

+ | (W - 7) 4, 2 ' 31 (£) - 18 4 3 - 4 ’ (C) 

(N - 7) e™ (() (<) +9 „!+ 2 > (<) <4+ 2 > (<) 

+18 < 4 W > (c) k 31 (<) - 9 k +‘>(<)| 2 


r (l,2) 


(C) 


+ (1 — s nm ) £ r 


18 C„ 


9 ( r , 


(N- 7)/2 + 2 <4 M) (C) 


(Cn - C m f (Cn - CJ^ 


—2 (3 N - 7) + (9/2) (IV - 7) <4 M) (C) - 9 (C) + 9 a ( n hl) (C) 


r (2,2) 


r(M) 


(C„ - C m ) 


(65c) 


and so on (via standard, but increasingly tedious, computations). 
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